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Abstract 

We show that the braid-group extension of the monodromy-based topolog- 
ical quantum computation scheme of Das Sarma et al. can be understood 
in terms of the universal R matrix for the Ising model giving similar results 
to those obtained by direct analytic continuation of multi-anyon Pfaffian 
wave functions. It is necessary, however, to take into account the projec- 
tion on spinor states with definite total parity which is responsible for the 
topological entanglement in the Pfaffian topological quantum computer. 



1 Introduction 

Quantum Computation [1] in general is expected to allow us to solve compu- 
tational problems that are hard to attack by classical methods. Unfortunately, 
this exponential speed-up over classical computation has not been demonstrated 
so far because of the overwhelming noise and decoherence due to the coupling 
between the qubits and their environment. The Topological Quantum Compu- 
tation (TQC) [2-4] is an intriguing proposal to use the braiding operations of 
non-Abelian quasiparticles in certain strongly correlated electron systems, such 
as the fractional quantum Hall (FQH) liquids, as quantum gates. The main ad- 
vantage of using this highly complicated approach is that the encoding of quan- 
tum information, as well as the representation of the corresponding quantum 
operations, are naturally immune against noise and decoherence due to the pres- 
ence of a bulk energy gap in the excitations spectrum. This has been called 
topological protection of quantum information and quantum operations because 
quantum information is naturally encoded in topological quantum numbers and 
then quantum gates must be implemented by topologically nontrivial operations. 
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One of the remarkable theoretical achievement in this direction was the TQC 
scheme of Das Sarma et al. [4] in which the qubit has been constructed by 4 
non-Abelian Ising anyons localized on 4 antidots in a Pfaffian FQH liquid and 
the qubit measurement was implemented electronic by Mach-Zehnder interfer- 
ometry. In this scheme the computational basis for quantum computation was 
realized in terms of 4-point correlations functions of the critical two-dimensional 
Ising model and a logical NOT gate has been implemented by a monodromy 
transformation. 

The monodromy-based TQC scheme of Ref. [4] has been extended [5, 6]to 
include braid-group operations by using the exchange matrices for 4-quasiholes 
Ising-model wave functions in the Pfaffian FQH state. Unfortunately it turns out 
that the representation of the (infinite) braid group B4 is actually finite and as a 
result not all single-qubit gates could be realized simply by braiding. In addi- 
tion to the single-qubit construction, an important two-qubit generalization has 
been proposed in Refs. [5, 6], in terms of 6 Ising anyons, and two-qubit gates 
have been implemented by braiding operations from the Ising-model represen- 
tation of the braid group Bq. The elementary braid matrices for 6 quasiholes 
have been computed explicitly from the 6-quasiholes Pfaffian wave functions 
and a completely topologically protected construction of the ControUed-Z and 
Controlled-NOT gates has been given [5, 6]. 

In this Report we will demonstrate that the above braid matrices could be 
generated from the universal R matrix for the Ising model, obtained before [7,8], 
by introducing an appropriate projector taking into account the fermion parity 
conservation in the Ising-model correlation functions. This projection turns out 
to be the origin of the topological entanglement which is characteristic for the 
Pfaffian topological quantum computer. 

2 Quantum computation in general 

The quantum bit of information, called the qubit, is a quantum state which be- 
longs to a two-dimensional complex vector space usually represented in an or- 
thonormal basis { 1 0) , 1 1) } as a normalized coherent superposition 

= a\0) + j3\l), where a, /3 e C and \af + \l3f = l. 

In any concrete quantum computation approach we have to first define the com- 
putational basis |0), |1) and then initialize the system by constructing e.g., 
\t{j) = |0) or = |1) (or some known superposition state). The measurement 
is understood as a standard orthogonal projectors family defining the projected 
state after measurement and the corresponding probability [1]. The multiple 
qubits are constructed as tensor products of the individual qubits so that, e.g., 
the n-qubit space of states is (the projective space of) 

The initialization of the multi-qubit system in quantum computation is usually 
done by preparing the system in a specific state in the multi-qubit Hilbert space 

010011000.. .01 ^ 1010011000.. .01) e W". 
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Any quantum computation could be performed by applying to the input multi- 
qubit state a finite sequence of quantum gates, which are represented by unitary 
operators over and the output of the computation which is obtained after 
measurement is classical data that could be written on paper. The exponential 
speed-up of quantum computation is due to the quantum entanglement and to a 
phenomenon called quantum parallelism [1]. 

Any quantum gate can be approximated [1], with arbitrary precision with 
respect to E{U, V) = max|^) 1 1 (i7 — 1 1/') 1 1 , by products of 3 universal gates: 
H, CNOT and T. These gates are called discrete set of universal quantum gates 
and can be written explicitly as follows 

1 0" 
10 
1' 
10 

where H and T are given in the basis { 1 0) , 1 1) } while CNOT is in the two-qubit 
computational basis {|00), |01), |10), |11)}. 

Unfortunately the decoherence and quantum noise due to local interactions 
destroying coherent phenomena and flipping |0) |1) create unavoidable ob- 
stacles to building an efficient scalable quantum computation platform. There 
are hopes that the quantum error-correcting algorithms could help creating ef- 
ficient quantum computers from qubits and quantum gates that are not perfect, 
however, these might not be a big help on the scale of 1000-qubit systems. That 
is why it has become increasingly important to look for systems which are in- 
trinsicaUy immune against noise and decoherence. The topological quantum 
computers are good candidates for almost noiseless processing of quantum in- 
formation. 

3 What is Topological Quantum Computation? 

Because the quantum noise and decoherence are presumably due to local interac- 
tions we could try to avoid them by encoding quantum information non-locally, 
for example, by using some topological characteristics of the system such as, 
homotopy classes of quasiparticles exchange paths. In this case quantum in- 
formation is inaccessible to local interactions, because they cannot distinguish 
between |0) and |1) and hence cannot lead to decoherence and noise, provid- 
ing in this way, a topological protection of qubit operations. Then the quantum 
gates should be some topologically non-trivial operations, such as changing the 
homotopy classes, which could be implemented by braid operations (exchanges 
of quasiparticles in the plane). 

However, in order to realize even a single topologically protected qubit, 
which belongs to a 2-dimensional space, we need state degeneracy in the plane 
and therefore higher dimensional representations of the braid group. The corre- 
sponding quasiparticles (known in the mathematical literature as plektons) are 
called non-AbeUan anyons and possess a number of very strange and interesting 
properties. 
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Perhaps the most promising non-Abehan candidate is the FQH state with 
filUng factor v = 5/2 that is routinely observed in ultrahigh-mobility samples 
[9, 10] in the second Landau level, and is believed to be in the universality class 
of the Moore-Read Pfaffian state whose two-dimensional conformal field theory 
X Ising contains the critical Ising model with central charge 1/2. 

The topological protection in this case comes from the fact that the residual 
noise and decoherence are due to thermally activated quasiparticle-quasihole 
(that might execute uncontrolled braidings), which are exponentially suppressed 
at low temperature (T 5 mK) by the bulk energy gap (A 500 mK) [9, 10]. 
In Ref . [4] the error rate has been estimated to be 

Error rate ~ ( exp (-j^) ~ lO'^^ 



which would lead to unprecedented precision of quantum information process- 
ing. 

4 Braiding: Extend TQC scheme of Das Sarma et al. 

The main idea of the TQC scheme of Ref. [4] is to use as a qubit the state de- 
scribed by the Pfaffian wave functions with 4 quasiholes. The quasiholes with 
coordinates r/i , . . . , 774 are localized on 4 antidots inside of an incompressible 
electron hquid formed by a large number N of electrons (or holes) with coordi- 
nates z-i,. . . ,zn and the corresponding wave function could be expressed as a 
chiral CFT correlation function 

N 

*4qh(?7l, Vi; Zi) =*(V'qh(»7l)V'qh(j72)V'qh(??3)V'qh(??4) V'h(^i)), (1) 

of the CFT operators representing the holes and quasiholes 

^^(z) = ^(z) : c'V^<?(^) : and V7qh(r/) = (7(7?) : e*^"^^"^ :, 

where and a{r]) are respectively the Majorana fermion and the chiral spin 
field of the Ising model and (j){z) is a normalized u(l) boson. The 4-quasihole 
wave function (1) have been computed explicitly in Ref. [7] to be 

where the functions ^'4°^^ <-> iO) , 1 1), which will serve as a computational 
basis for the Pfaffian qubit, have the form {riab = Va — Vb) 



*4qh (^l''72,'73,?74;2:i, ■■■,Zn) = , i*(13)(24) ± W(i4)(23) j 



and are expressed in terms of the single-valued Pfaffian functions [6, 7] 

,T, pr f {Zi - Va) {Zi - Vb) {Zj " Vc) {Zj - Vd) + {i j) 



Yl {zt - z/f, {a<b, c<d) 

l<i<j<N 
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and X = ?7i4r/23/(T/i3?724) is the anhamonic CFT invariant crossratio. 

The fermion parity is a subtle issue in the Ising model [6, 1 1]. The Ramond 
superselection sector, is double degenerate as a result of the presence of the 
Majorana fermion zero mode, i.e., there are two chiral primary fields a± with 
CFT dimension 1/16, whose subscript is the chiral fermion parity 

[■^0,7f]+ = 0, V'o = ^> 7f = 1 7FCr±7F = ±£7±. 

The physical quasihole, however, is represented by a{r]) = {(T+{r]) + 
(*?))/ because of the modular invariance requirement and the necessity of 
the so called GSO projection [11], which is actually the origin of its non-Abelian 
statistics expressed in the fusion rule a x a = 1 + ^. 

The conservation of the fermion parity implies that e.g., the correlation func- 
tions involving 4 a fields and (possibly) an even number of Majorana fermions 
are non-zero (CTeiCre^creafTei) 7^ only if 6162 = 6364. On the other hand it 
follows from the Ising-model operator product expansion that cr+(j+ ~ (t_(j_ 
and (7-1- tT_ ~ tT-fT-i-, which together with 6162 = 6364 reduces the dimension of 
the 4cr-correlation functions space to dim H4cr = 2. This allows us to write the 
computational basis in the CFT description as 

|0) = {a+a+c7+a+), |1) = {a+a-a+a-), 

i.e., the state of the qubit is determined e.g., by the second pair fT,.;, o",;., of a fields 
and the states of other pair is fixed by fermion parity conservation. In general [7] 
the correlation functions space dimension for 2n a fields is dim 'H2na = 2"~^, 
which means that 2n Ising anyons could be used to represent n — 1 qubits. 

5 Braid matrices from tlie Pfaffian wave functions witli 4 quasi- 
lioles 

Now that we know how to write the computational basis in terms of Ising model 

correlation functions we could represent quasiparticle braidings, that would be 
used as quantum gates, as analytic continuation of their coordinates. For exam- 
ple the counterclockwise braiding of the quasiparticles with coordinates rji and 
ri2 could be executed by the analytic continuation along the circle by 

> m+Jl^^i^tm^ ,^rn+^_^,,,^^ ^^^^^ 
/I 2 2^2 2 

Doing carefully the analytic continuation of the Ising-model 4-point functions 
-\/l ± \/x, with X = 77147/23 / (7/i3'724), and taking into account the contribution 
of the single-valued functions ^'(a6,cd) we could obtain [6, 7] the generators of 
the (positive parity) 2-dimensional representation of the braid group B4 to be 



p(4) _ d(4) 
-"■12 — ^34 



1 
i 



•-23 



V2 



1 -i 
-i 1 



(2) 



Alternatively, these generators have been obtained in Ref. [12] by using only 
quantum-group methods, which is more compact and perhaps more elegant. 



6 



Topological Quantum Computation 



However, the direct derivation [6] from the multi-quasihole Pfaffian wave func- 
tions is necessary for many reasons. 

It turns out that the representation generated by Eq. (2) is of finite order [6]. 
Given the explicit form of the generators of a finite group, we could enumerate 
all group elements by using an optimized Dimino algorithm [13]. For the repre- 
sentation of the group B4 generated by (2) the Dimino's algorithm gives (using 
Maple 10 with modules) 

[Image (64) | = 96. 

For general 2n Ising anyons the representation of the braid group B2n (with 
2n > 6) is finite as well and its order has been given in Ref. [14] 

IT M f 22"-i(2n)! for n = even 

|Image(S2n)| = | 22n(2n)! for n = odd " 

This is too bad because for universal TQC we need Image {B^n) to be dense in 
the unitary group. This means that not all quantum gates could be realized by 
braiding, hence, not all of them would be fully topologically protected. Nev- 
ertheless, most of the important quantum gates, especially those from the Clif- 
ford group, can be implemented by braiding [6]. For example the single-qubit 
Hadamard gate can be constructed [6] by 3 elementary exchanges as shown in 
Fig. 1. 

[1 1" 

V2 V- -1. 




Figure 1. Hadamard gate implemented by braiding 4 Ising anyons. 



6 Exchange matrices for 6 Ising quasiholes: tlie two-qubit gates 

The two-qubit states could be represented by Pfaffian wave functions with 6 
quasiholes localized on antidots [6] and the computational basis would be e.g., 

|00) = ((T_|_cr_|_cr_|_(T_|_cr_|_(T_|_), |01) = (cr_|_(T_|_cr_|_(T_cr_|_cr_) 
1 10) = {cr+a-a+a-a+(j+) , |11) = (cr_|_cr_<T+cr+cr+cr_). 

In order to obtain the braid matrices for exchanging 6 Ising anyons, which 
would represent the two-qubit gates, we don't need to consider the full 6- 
quasiholes Pfaffian wave functions which are not even known explicitly. Instead, 
we could use the general expression similar to Eq. (1), however with 6 quasihole 
operators V'qh inside the correlator. Then we can use the fact that the braiding 
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operations are independent of the distance between the particles, especially for 
those which are not involved in the exchange, and we could fuse them to ob- 
tain a 4 quasihole Pfaffian wave function for which we could apply the results 
in Eq. (2). For example, if we want to compute R^f!] then we could first fuse 
% in which case the two-qubit computational basis could be written as 

|00) {a+a+a+a+), |01) {a+a+a+a-^p) 



1 10) {a+a-a+a-), 

'75— >'76 



1 11) ^ {a+a-a+a+ip). 

'75— >'76 

?(4) 



Then the exchange r]i ^ r]2 is represented by and it is easy to see that 

Similarly, by first fusing rj^ ^ rje, we could compute R^^ as follows 





R'. 



(6) 
■23 



C 4 

71 



1 -i 

1 

-i 1 

-i 



R^^ O I2 



In the same way, however, by fusing first rji 



R 



(6) 
45 



e'4 

7f 
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—i 








—i 
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—i 
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1 



r]2, we could obtain 

= I2 0i^it^ 



and 



=diag(l,2,l,?:) = l2®i?f^ 



34 



Unfortunately, this approach doesn't work for R^^^ because, as we will see later, 
it is not factorizable due to the fact that the states of the quasiholes at r]s and 774 
depend on both qubit's states. However, due to a general argument [6] based on 
the superselection rule in the Neveu-Schwarz sector, R^J must be diagonal [6]. 
Therefore, we could obtain R^^ directly from the OPE for 773 774 

|01) vli^if^+c^+^f^+c^-) 



|00) »734^^^(o'+o-+o'+o'+)> 

'73— ►'74 



|10) 



'73— >'74 



%f (0-+0--V'0-+CT+), |11) 



»73— >»74 



-1/8, \ 
%4 (0-+f^-0-+0--) 



Taking into account the extra factor rj^^^ coming from the AbeUan part of the 
quasiholes, and executing the exchange by 7734 e^'^iju we obtain the exchange 
matrix for 773 ^ 774 to be 



-R^4^ = diag(l,i,i, 1). 



(3) 



It is worth stressing that the 6-quasihole exchange matrix r'^^ cannot be ex- 
pressed as a factorized tensor products of the 4-quasiholes exchange matrices 
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and the unit matrix I2. This important fact, which is responsible for creating 
entanglement purely by braiding, will be explained later in the universal R ma- 
trix approach. It is this topological entanglement which allows us to implement 
entangling two-qubit gates, such as the CNOT gate, purely by braiding 6 Ising 
anyons. Indeed, using the explicit formulae for the 6-quasiholes braid matrices, 
it is easy to check that the CNOT can be represented by 7 elementary exchanges 
as follows [5, 6] (see Fig. 2 for the corresponding braid diagram) 

CNOT = R^l Rn^R-iiRi2Rb6RibRsi ■ 




CO 



-0- 



Figure 2. Braid diagram for realization of the CNOT gate by braiding 6 Ising anyons. 
The symbol to the right is the standard notation for CNOT. 



7 The universal R matrix approach 



The quantum group structure for the Pfaffian model has been analyzed in 
Ref. [12]. Using the affine coset description, Ising= sm(2)2/w(1), one could 
identify its quantum group as Uq{sl{2)) with q = e~'^. 

As we have already mentioned the n-qubit states are represented by 2(n + 1) 
Ising anyons which reaUze [7] spinor irreducible representations of SO{2n + 2) 
of dimension 2". Therefore the spinor representation of the braid-group gener- 
ators could be expressed in terms of the spinor representation of the universal 
i?-matrix for Uq{sl{2)) with q = e~'7 as follows 



Rk,kt? = 12 ® • • • ® I2 O-R O I2 O ' 



(4) 



k-l 



The spinor representation of the R matrix for the Pfaffian FQH state (or for the 
Ising model) depends on the parity of k [7, 8, 12]. When k is odd the R matrix 
acts only on the A;-th factor and has the form 

"1 0" 



R- 







while for even kthe R matrix acts on the fc-th and {k + l)-th factor as 



R 



1 



1 
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Notice however, that the exchange matrix (4) has dimensions 2^"+^ x 2^"+^ 
while the spinor representations of SO(2n + 2) have dimension 2". In order 
to obtain the action of the (irreducible) exchange matrices in the spinor repre- 
sentation with positive parity we have to apply the corresponding projector For 
example, in the case of 6 Ising anyons, using the obvious notation <j+cr+ ~ + 
and cr+f7_ ~ — , we have to project the states 



! + ++),! + +-), ! + -+),! + — ), I - ++), I - +-), I - -+), I - — ) 



to the two-qubit basis { I + ++) , |H ), | — H — ), | h) } with positive parity. 

Therefore, to obtain the braid-group generators in the representation we have to 
apply the projector 



Now let us go back to the braid matrix for the exchange of 773 and 774 in the 6- 
anyon representation. Applying the projector P2 and deleting all nuU rows and 
columns we get the same result as Eq. (3) 



Obviously factorizability of the braid matrix (inside the parentheses above) is 
lost after projection and this projection is the origin of the topological entangle- 
ment discussed in Sect. 6. 

It is worth mentioning that the Artin's relation for the spinor representations 
of the braid-group generators follow from the Yang-Baxter equation for the R 
matrix only if the projectors commutes with I2 • • • 12 -R 12 • • • 12, 
however, it turns out that for the Ising anyons this is always the case. 

Finally it might be interesting to note that this representation of the braid 
matrices in terms of the R matrix is a nice illustration of the idea of Kauff- 
man et al. to use universal R matrices for TQC [15]. Of course we have to 
be aware that some specific details, such as parity projectors creating topolog- 
ical entanglement and the impossibiUty to reaUze by braiding all single-qubit 
gates, may be model-dependent, however, we have to admit that the opportunity 
for topological quantum entanglement due to braiding of non-Abehan anyons is 
remarkable. 

8 Conclusions 

Several open problems should be mentioned: first of all it might be very useful 
if we could construct the missing 7r/8 gate by braiding plus something else in 
such a way to obtain maximal protection from noise. Alternatively one may try 
to construct a protected Toffoli gate. The embedding of all single-qubit and two- 
qubit gates into three-qubit system still has to be resolved. Another question is 
how to approximate the quantum Fourier transform more efficiently and finally 
it would be necessary to analyze the possible error sources in this TQC scheme. 



P2 :=diag(l, 0,0, 1,0.1,1.0). 
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